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ABSTRACT 

The  nucleate  boiling  phenomena  is  analysed  for  Ti02  and  Al203  aqueous  Nanofluids  using  computational  fluid 
dynamics.  The  use  Nanofluid  as  a  working  fluid  significantly  enhances  the  boiling  critical  heat  flux  (CHF). 
It  has  been  found  that  the  CHF  enhancement  in  boiling  is  dependent  on  the  type  of  Nanofluid  as  well  as  its 
concentration.  The  present  study  observed  that  by  increasing  the  nanoparticles  concentration,  the  bubble  frequency 
increases.  Two  sets  of  Nanofluids,  viz.  Al203  and  Ti02  with  Volume  of  Fraction  (VOF)  of  0.05,  0.1,  and  0.15  have  been 
considered  for  the  study.  The  obtained  simulation  results  show  that  by  increasing  the  nanoparticles  concentration,  the 
TSHF  and  HTC  increase  proportionally.  Ti02  nanoparticles  with  various  VOF  give  better  results  after  a  time  interval  of 
1.4  s  compared  to  Al203  and  water. 
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NOMENCLATURE 

(p  =  Particles  volume  concentration 
Pf  =  Density  of  the  base  fluid,  kg/m3 

р, p=  Density  of  the  Nanoparticles,  kg/m3 

Nanofluid  viscosity,  kg/ms 
Water  viscosity,  kg/ms 

с, jrf  =  Heat  capacity  of  the  Nanofluid,  J/kg°C 
Cyjj  *=  Heat  capacity  of  the  base  fluid,  J/kg°C 

CpfJ=  Heat  capacity  of  the  Nanoparticles,  J/kg°C 

kp  =  Thermal  conductivity  of  the  Nanoparticle,  W/m°C 
kbj=  Thermal  conductivity  of  the  water,  W/m°C 
<pv=  Volume  fraction 
Pe  =  Peclet  number 
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Re  =  Reynolds  number 

Pr  =  Prandtl  number 

ABBREVIATIONS 

VOF  =  Volume  of  Fraction 

HTC  =  Heat  Transfer  Coefficient 

TSHF  =  Total  Surface  Heat  Flux 

INTRODUCTION 

Boiling  heat  transfer  plays  an  important  role  in  a  variety  of  industrial  process  and  applications  like  refrigeration, 
heat  exchangers,  cooling  of  electronic  components  and  power  generation.  Boiling  heat  transfer  enhancement  process  is 
vital  and  could  make  a  significant  impact  on  the  industrial  growth. 

Nanofluids  are  one  of  the  new  kind  of  heat  Transfer  fluids  in  which  Nano  meter  sized  particles  are  dispersed 
uniformly.  Recently  many  researchers  have  proposed  a  new  approach  that  will  considerably  enhance  the  boiling 
phenomena.  In  the  study  of  boiling  regimes,  nucleate  boiling  is  an  efficient  heat  transfer  mechanism.  However,  critical  heat 
flux  (CHF)  phenomenon  is  the  thermal  limit  during  the  heat  transfer  phase  change  where  heat  transfer  is  maximized 
followed  by  drastic  degradation.  The  study  of  bubble  nucleation,  growth  and  detachment  are  inherently  complex 
phenomena,  characterized  by  formation  of  vapour  bubbles  on  heated  wall.  Numerical  simulation  for  analysing  the  bubble 
dynamics  and  heat  transfer  in  nucleate  boiling  has  been  studied  by  several  investigators  [1-7].  Recently,  numerical  studies 
have  been  carried  out  for  analysing  bubble  dynamics  and  heat  transfer  on  a  horizontal  surface.  Use  of  Nanofluids  play  a 
key  role  in  improving  the  process  efficiency,  increasing  the  CHF  and  reducing  the  operating  cost;  this  is  where  Nanofluids 
emerged  as  efficient  heat  transfer  fluids.  Research  on  enhancement  of  CHF  using  Nanofluid  under  convective  flow 
condition  has  also  been  studied,  but  to  a  lesser  extent.  Interestingly  majority  of  the  experimental  data  are  on  the  effect  of 
Nanofluid  on  the  CHF.  There  is  a  significant  gap  in  the  study  of  bubble  nucleation  and  detachment  which  also  plays  a  vital 
role  in  studying  the  heat  transfer  applications.  The  main  advantage  of  using  Nanofluid  for  the  heat  transfer  applications  is 
the  ability  of  the  Nanofluids  to  alter  these  properties  like  thermal  conductivity,  specific  heat,  density,  surface  wettability. 
These  properties  can  be  adjusted  by  varying  the  Nanoparticles  concentrations  in  the  base  fluid. 

In  the  present  work,  two  metal  oxides  -  Nanoparticles  A IX)  and  Ti02  are  considered  for  the  analysis,  which  are 
suspended  in  water.  Numerical  simulation  has  been  carried  out  for  analysing  the  bubble  dynamics,  bubble  detaching 
frequency,  surface  critical  heat  flux  and  heat  transfer  coefficient  and  the  results  are  compared  with  water  under  transient 
condition.  Kim  and  Kim  (2009)  [1]  reported  that  Nanoparticle  motion  depends  on  hydrodynamic  forces  formed  around  the 
surface,  which  lead  to  disappearance  of  influence  of  Nanoparticle  and  completely  depend  upon  the  volume  concentration 
(vol  %).  They  conducted  experiments  on  nucleate  boiling  using  different  TiCB,  AUO3  Nanofluids  and  stated  the 
enhancement  of  CHF  up  to  100%  when  compared  with  water.  They  stated  the  reason  as  deposition  of  oxides  of 
Nanoparticles  coated  on  heater  surface,  which  improves  the  wettability,  which  in  turn  enhances  CHF.  Experimental 
investigations  were  conducted  by  You  et  al.  [2]  on  a  flat  square  plate  heater  placed  in  Nanofluid  with  a  fine  proportion  of 
Nanoparticles  mixed  and  observed  enhancement  of  CHF. 

Li  (2016)  [3]  conducted  CFD  analysis,  observed  velocity  vector  contours  and  stated  that  the  vapour  velocity 
during  the  departure  is  much  more  than  the  vapour  velocity  during  the  growth  of  bubble.  He  has  shown  the  Eddies 
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formation  and  flow  pattern  of  liquid  around  the  bubble  and  inside  the  bubble.  Buongiorno  and  Hu  (2010)  [4]  after 
conducting  experimental  investigations  have  observed  the  Enhancement  of  CHF  by  Nanofluid  where  nucleate  boiling  takes 
place.  They  also  stated  the  incremental  CHF  for  both  A1203  and  SiC>2  of  200%  and  60%  respectively  for  a  small  proportion 
Nano  particle  dispersed  in  water  and  compared  with  water. 

Wei  et  al.  (2015)  [5]  found  that  the  effect  of  growing  diameter  is  inversely  proportional  to  gravity.  They  came  to 
the  conclusion  using  the  thermal  lattice-Boltzmann  method  and  correlation  between  the  departure  diameter  and  the  gravity. 
The  bubble  diameter  at  the  time  of  departure  is  given  by: 

£7 

y(Pi  - 

After  conducting  both  experimental  and  CFD  analysis  on  pool  boiling  by  taking  a  heated  wire,  using  Nanofluid  as 
fluid  domain  Girish  Sapre  et  al.  (2018)  [6]  concluded  Ti02  has  a  high  bubble  frequency  compared  with  A1203  which  result 
increase  in  HTC  and  TSHF.  After  theoretical  study  author  Narendra  et  al.  (2018)  [7]  has  observed  that  Nanofluid  thermal 
conductivity  increases  with  increase  in  volume  concentration  of  nanoparticles,  comparing  the  thermal  conductivity  of  both 
Ti02  and  A1203  Nanofluid  concluded  Ti02has  high  thermal  conductivity  then  A1203.  ANSYS  Fluent  software  package  has 
been  utilized  to  solve  many  thermal  engineering  problems  successfully  and  many  of  the  results  are  found  to  be  reasonably 
in  good  agreement  with  test  data  [8-11].  Base  fluid  as  well  as  NanoFluids  characteristics  have  been  examined  in  [12-16]. 

THERMO  PHYSICAL  PROPERTIES  OF  NANOFLUID 

Akilu  et  al.  (2016)  [17]  has  given  the  correlation  for  obtaining  thermo-physical  properties  for  Nanofluids  after 
considering  all  the  parameters  like  VOF,  particle  size  and  shape  of  the  nano  particles.  Koca  et  al.  (2018)  [18]  has  stated  the 
effect  of  nano  particle  size  on  the  viscosity  of  Nanofluid,  which  will  affect  the  bubble  growth  and  bubble  frequency.  The 
following  equations  are  used  to  calculate  the  properties  of  Nanofluids. 


Pnf  =it-  p>p bf + <Pp p 
Where; 

=  Density  of  the  Nanofluid,  kg/m3 

ip  =  Particles  volume  concentration 
_  Density  of  the  base  fluid,  kg/m3 
py  -  Density  of  the  Nanoparticles,  kg/m3 

•  Pnj  =  I  2'E-<p')'rtbf 
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Where; 

flnf=  Nanofluid  viscosity,  kg/ms 
Water  viscosity,  kg/ms 

Cpnf  pnf 

Where; 

Cp.nf  =  Heat  capacity  of  the  Nanofluid,  J/kg°C 
£phf=  Heat  capacity  of  the  base  fluid,  J/kg°C 
C  =  Heat  capacity  of  the  Nanoparticles,  J/kg°C 

jr r  _  , 

tep+2kj/--Ckp-ki/>pp  bf 

Where; 

kp=  Thermal  conductivity  of  the  Nanoparticle,  W/m°C 
kbf=  Thermal  conductivity  of  the  water,  W/m°C 
%=  Volume  fraction 

•  Nu  =  0.4328  (l+11.285^)(1'7S4'Fe(I,-2iS)i;ea'333Pr(1-4 
Where, 

<p  =  Volume  fraction 
Pe  =  Peclet  number 
Re  =  Reynolds  number 
Pr  =  Prandtl  number 

Following  tables  [1-3]  give  the  calculated  property  values  of  different  Nanofluids: 


Table  1:  Properties  of  Al203-Water  Nanofluid 


A1203- Water 

Volume  Fraction 

Density  [kg/m3] 

Cp  [J/kg-K] 

p  [kg/m-s] 

K  [W/m-K] 

0.05 

1145.745 

3589.25 

9.585E-4 

0.70514 

0.1 

1294.39 

3134.965 

1.065E-3 

0.807006 

0.15 

1443.035 

2774.264 

1.1715E-3 

0.9202 
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Table  2:  Properties  of  Ti02-Water  Nanofluid 


Ti02- Water 

Volume  Fraction 

Density  [kg/m3] 

Cp  [J/kg-K] 

p  [kg/m-s] 

K  [W/m-K] 

0.05 

1147.24 

3574.421 

9.525E-4 

0.6904 

0.1 

1297.39 

3109.776 

1.065E-3 

0.7748 

0.15 

1447.535 

3036.229 

1.1715E-3 

0.86696 

Table  3:  Properties  of  Water 


Water 

Density  [kg/m3] 

Cp  [J/kg-K] 

p  [kg/m-s] 

K  [W/m-K] 

998.2 

4182 

0.6 

0.001003 

COMPUTATIONAL  DOMAIN 


A  rectangular  section  of  0.004m*0.008m  computational  domain  is  set  up  for  bubble  capture,  resembling  a 
cylindrical  vessel  taken  in  2D. 

MESH  METHOD 


Grid  Independence  studies  have  been  conducted  to  acquire  consistent,  accurate  and  reliable  results.  It  is  observed 
that  the  number  of  elements  and  nodes  kept  on  increasing  until  the  heat  transfer  coefficient  (HTC)  of  water  became  stable. 
The  summary  of  numerical  experiments  conducted  for  deciding  the  grid  size  has  been  tabulated  below  (Table  5).  A  graph 
is  also  plotted  to  show  the  grid  independence. 


Table  4:  Results  of  Grid  Independence  Studies 


No.  of  Elements 

Heat  Transfer  Coefficient 
of  Water  at  0.2  s 

14323 

1.842 

16201 

1.876 

18822 

1.985 

20301 

2.030 

20811 

2.029 

30242 

2.031 

Elements 


Figure  1 

Based  on  the  above  numerical  experiments  conducted  and  to  maintain  less  computation  time  and  space,  the 
number  of  elements  has  been  fixed  at  20301  and  the  number  of  nodes  at  20,000.  The  mesh  has  been  generated  by  using 
ICEM  CFD  (meshing  tool). 
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Figure  2  represents  the  mesh  and  boundary  conductions  and  exponential  grid  spacing. 


condenser  wall  273K 


Heater  400K 


Figure  2:  Exponential  Gird  with  Representation  of  Boundary  Conductions 


Micro  level  Spacing  of  0.0001  mm  is  given  at  the  heater  surface  to  capture  the  bubble  dynamics  and  bubble 
departure.  Bi-exponential  grid  ratio  is  given  as  1 .5  as  shown  in  Figure  1 . 


GOVERNING  EQUATIONS 


The  analysis  is  based  on  transient,  two-phase,  three-dimensional,  turbulent  flow  with  the  Boussinesq 
approximation  for  the  buoyancy  term  and  VOF  model  for  capturing  the  bubble  dynamics  under  multi -phase.  The  governing 
equations  run  in  the  background  of  the  software  interface;  these  equations  describe  continuity,  momentum  and  energy 
transport.  The  following  governing  equations  are  used  in  computation  fluid  dynamics  Xiaobin  et  al.  (2015)  [19]. 

Continuity  Equation 

£(°,p,)+»<w.)  (1) 

P  =  Liquid  phase 
Q  =  V  apour  phase 


a-  =  Volume  fraction 


vq=  Velocity  of  phase  q  (vapour) 

Momentum  Equation 

+  V  ( ocqPq vTtr*)  =  —acVP  +  V  r  +  ocqPq^+  Rpq  +  (  Flift  +  Fcd  q  ) 

Energy  Equation 


(2) 


d_ 

Dt 


0«/lA  )  ■+■ v  0«/vcA ) 


=  oc. 


dt 


-F  rV  uq  —  V  cfq 


P= 1 


(3) 
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Bottom  edge  as  heater  is  provided  with  a  constant  temperature  of  400  K,  top  edge  as  condenser  wall  is  provided 
with  a  constant  temperature  of  273  K,  remaining  edges  are  set  adiabatic  conditions. 

RESULTS  AND  DISCUSSION 


Initial  bubble  growth  and  departure  time  are  numerically  investigated  in  this  investigation.  Horizontal  heating  is 
applied  to  the  bottom  wall  surface  by  maintaining  a  constant  temperature  of  400  K.  Similarly,  top  wall  is  kept  at  a  constant 
temperature  of  273  K  (Condensation  condition)  and  the  other  walls  are  given  adiabatic  conditions. 


To  observe  the  bubble  growth  and  bubble  dynamics,  small  time  step  of  lxlO"1  s  is  taken  during  numerical 
simulation.  Hence,  simulated  initial  bubble  along  with  time  is  observed  from  image  obtained  from  Numerical  study  for 
water,  Al203-Water  Nanofluid  and  Ti02-Water  Nanofluid  at  three  different  VOFs  of  0.05,  0.1  and  0.15  [Figures  2 -4c]. 
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Figure  2:  Phase  Contour  for  Water  During  Bubble  Departure  (Red  Colour 
Represents  the  Liquid  Phase  and  the  Blue  Represents  the  Vapour  Phase) 

In  Figure  2  indicates  the  bubble  growth  and  detachment  in  a  heated  surface  simulation  of  the  bubble  detachment 
has  been  carried  out  at  the  centre  of  the  compaction  demine  as  time  evolves  the  bubbles  gradually  groups  from  the  heated 
surface  and  it  will  moves  upwards  due  to  gravity  and  the  density  difference  there  by  bubble  experiences  an  upward  force 
leading  to  its  shape  deformation  and  lift  off  takes  place  from  the  wall  surface  Xu  et  al.  (2018)  [20].  As  the  bubble  moves 
away  from  the  wall  its  shape  deforms  again  due  the  unbalance  action  caused  by  bouncy  and  viscous  drag  force  from  the 
surrounding  liquid. 
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Figure  3a:  Phase  Contour  for  A1203- Water  Nanofluid  with  a  VOF  of  0.05  During 
Bubble  Departure  (Red  Colour  Represents  the  Liquid  Phase  and  the 
Blue  Represents  the  Vapour  Phase) 


■ 

U 


Contours  of  Volume  fraction  (al2o3_0.1)  (Time=1 ,4300e+00)  Mar  03.  2019 

ANSYS  Fluent  15.0  (2d.  dp,  pbns.  vof.  ske,  transient) 


Figure  3b:  Phase  Contour  for  Al203-Water  Nanofluid  with  a  VOF  of  0.1  During  Bubble  Departure 
(Red  Colour  Represents  the  Liquid  Phase  and  the  Blue  Represents  the  Vapour  Phase) 
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Contours  of  Volume  fraction  (al2o3-0. 15)  (Time=1.7800e+00)  Mar  01,  2019 

ANSYS  Fluent  15.0  (2d.  dp.  pbns.  vof.  ske.  transient) 


Figure  3c:  Phase  Contour  for  Al203-Water  Nanofluid  with  a  VOF  of  0.15  during  Bubble  Departure 
(Red  Colour  Represents  the  Liquid  Phase  and  the  Blue  Represents  the  Vapour  Phase) 


Figure  3a  indicates  the  bubble  formation  and  departure  for  the  A1203  Nanofluid  with  0.05  VOF  and  Figures  3b 
and  3c  indicate  bubbles  formation  and  departure  for  the  same  Nanofluid  and  with  0.1  and  0.15  VOF.  In  all  the  three  figures 
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we  can  observe  that  by  increasing  the  percentage  of  Nanoparticles  heat  flux  increases  and  bubble  departures  from  surface 
with  a  time  gap  and  the  waiting  time  for  the  formation  of  the  new  bubble  decreases,  it  means  by  increasing  volume 
concentration  the  bubbles  form  rapidly  from  the  heated  surface  [21-22].  By  increasing  volume  concentration  the  bubble  lift 
of  diameter  decreases  steeply  due  to  increase  of  the  heat  flux. 


Figure  4a:  Phase  Contour  for  Ti02-Water  Nanofluid  with  a  VOF  of  0.05  during  Bubble  Departure 
(Red  Colour  Represents  the  Liquid  Phase  and  the  Blue  Represents  the  Vapour  Phase) 
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Figure  4b:  Phase  Contour  for  Ti02- Water  Nanofluid  with  a  VOF  of  0.1  during  Bubble  Departure 
(Red  Colour  Represents  the  Liquid  Phase  and  the  Blue  Represents  the  Vapour  Phase) 
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Figure  4c:  Phase  Contour  for  Ti02-Water  Nanofluid  with  a  VOF  of  0.15  during  Bubble  Departure 
(Red  Colour  Represents  the  Liquid  Phase  and  the  Blue  Represents  the  Vapour  Phase) 
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Figures  4a,  4b,  4c  indicate  the  bubble  formation  and  departure  from  the  heated  surface  for  the  Ti02  Nanofluid  for 
0.05,  0.1,  0.15  volume  concentration.  While  comparing  the  thermal  properties  of  A1203  and  Ti02  mono  fluids  it  has  been 
observed  that  the  thermal  conductivity  is  higher  for  the  A1203  Nanofluids,  therefore  we  can  observe  a  good  number  of 
bubbles  in  A1203  Nanofluids  Li  et  al.  (2016)  [23].  Whereas  the  bubble  frequency  is  more  for  tio2  it  means  the  consecutive 
bubble  formation  will  be  more  in  the  Ti02  Nanoparticles  due  to  the  lesser  specific  heat  capacity  compared  to  the  A1203 
Nanofluids.  The  bubble  departure  is  a  function  of  the  gravity  force  which  is  proportional  to  the  g0'5;  it  means  initial  bubble 
radius  influences  the  path  of  the  early  growth  period  which  is  going  to  reflect  in  the  simulation  for  a  great  departure.  Figure 
5  shows  bubble  dynamics  and  growth  of  bubble  and  detachment  of  bubble  from  the  heated  surface.  Time  steps  are  also 
provided. 


Table  5 


0.9(s) 

1.1(8) 

1.2(8) 

1.3(8) 

1.4(8) 

1.59(s) 

1.7(8) 

1.8(8) 

Figure  5:  Phase  Contour  of  Development  of  Bubble  and  Departure  of 
Bubble  form  Surface  for  A1203  of  VOF  of  0.1 


The  following  bubble  growth  and  departure  shown  in  the  figure  is  compared  with  the  results  obtained  by  the 
author  Bang  et  al.  (2005)  [24],  the  departure  time  steps  have  been  considered. 

Nanofluids  with  0.05%  of  alumina  show  similar  TSHF  when  compared  to  0.1%  and  0.15%  of  concentration,  due 
to  the  deposition  of  the  Nano  particle  on  the  heater  surface,  which  affects  the  boiling  phenomena  (CHU  et  al.,  2011)  [25]. 
Nanofluids  used  in  the  study  have  contributed  to  the  enhancement  of  TSHF.  Various  researchers  criticized  that  deposition 
of  Nanoparticles  on  the  surfaces  is  the  main  cause  beyond  the  TSHF  enhancement  due  to  the  effect  of  surface  liquid 
wettability  and  morphology  [26-27].  The  phenomenon  of  wettability  varies  depending  on  the  Nanoparticles  concentration 
with  the  base  fluid.  In  the  pool  boiling  phenomena,  TSHF  varies  with  various  concentration  of  Nano  particles.  At  a 
concentration  of  0.05,  the  critical  heat  flux  for  all  the  Nanofluids  increase  steeply  (Suriyawong  et  al.,  2013)  [27].  By 
increasing  the  Nanoparticles  concentration  it  has  been  observed  that  there  is  no  sizable  increase  in  the  TSHF.  All  the 
Nanofluids  seems  to  have  the  same  characteristics  of  TSHF  enhancement  qualitatively  (Ivey,  1967  and  Kim  et  al.,  2007) 
[28-29]. 

Figures  6  and  7  give  the  comparison  of  TSHF  Vs  HTC  for  A1203  and  Ti02  concentrations  at  different  time 

intervals. 
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Figure  6:  Comparison  of  HTC  Vs.  TSHF  of  A12Q3  (Vol%  of  0.05, 0.1, 0.15) 


Figure  7:  Comparison  of  HTC  Vs.  TSHF  of  Ti02  (Vol%  of  0.05, 0.1, 0.15) 

As  could  be  seen,  there  is  a  linear  variation  for  all  the  A1203  and  Ti02  concentration.  Suriyawong  et  al.  (2013) 
[30]  has  also  obtained  similar  results  during  their  investigation  with  Al203-water  Nanofluid. 

Figure  8  shows  fluctuation  of  HTC  along  with  time  on  the  heated  surface  when  investigated  and  compared 
between  Water,  Al203-0.05,  and  TiO2-0.05  Nanofluid. 


Figure  8:  Comparison  of  HTC  Vs.  Time  of  Water,  Al203-0.05,  TiO2-0.05 

Figure  8  shows  variation  of  HTC  with  respect  to  the  bubble  dynamics  time.  The  graph  shows  that  the  HTC 
decreases  between  0. 2-0.4  s,  and  remains  constant  from  0.4  s  onwards  for  Ti02  Nanofluid  with  0.05  concentration.  As 
could  be  seen  from  the  figure,  higher  heat  transfer  coefficients  have  been  observed. 
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Due  to  the  lesser  cp  value  of  3574.421  [J/kg-K],  when  compared  to  A1203  at  0.05  concentration  3589.25  [J/kg-K] 
from  the  property  it  has  been  observed  that  density  of  Ti02  Nanofluid  with  0.1  %  concentration  has  a  lesser  value 
compared  to  AI2O3  due  to  this  the  bubble  dynamics  for  the  Ti02  Nanofluids  is  efficient. 

Figure  9  shows  fluctuation  of  HTC  along  with  time  on  the  heated  surface  when  investigated  and  compared 
between  Water,  Al2O3-0.1  and  TiO2-0.1  Nanofluids. 


Figure  9:  Comparison  of  HTC  Vs.  Time  of  Water,  Al2O3-0.1,  TiO2-0.1 

Figure  9  shows  the  HTC  vs  time  for  A1203  and  Ti02  Nanofluids  for  0.1%  of  concentration,  compared  with  the 
water.  As  can  be  observed  that  form  a  time  period  of  0.2  to  0.4  s,  there  is  a  sudden  decrease  in  the  HTC  value  and  from  0.2 
to  0.4  s  and  the  HTC  will  remain  nearly  constant  from  0.6  s  onwards.  This  variation  is  captured  in  the  contour  plots  and  it 
has  been  observed  that  during  0.2-0. 4  s,  the  convection  current  plays  a  significant  role  and  the  small  bubbles  raises  from 
the  surface  and  collapse  due  to  the  forces  of  buoyancy  from  the  time  period  0.6s  onwards.  It  can  also  be  seen  that  bubbles 
rise  constantly  throughout  the  structure  in  a  uniform  manner;  therefore  heat  transfer  coefficient  will  be  nearly  constant 
(Kamatchi  et  al„  2015)  [31]. 

Figure  10  shows  fluctuation  of  HTC  along  with  time  on  the  heated  surface  when  investigated  and  compared 
between  Water,  AI2O3-O.I5  and  TiO2-0.15  Nanofluids. 


Figure  10:  Comparison  of  HTC  Vs.  Time  of  Water,  AI2O3-0.15,  TiO2-0.15 


Impact  Factor  (JCC):  7.6197 


SCOPUS  Indexed  Journal 


NAAS  Rating:  3.11 


The  Numerical  Analysis  of  Bubble  Growth  and  Bubble 
Frequency  in  Nucleate  Boiling  Using  Nanofluids 


1271 


Figure  10  shows  that  for  a  time  period  up  to  0.4  s  the  HTC  decreases  steeply  for  all  the  three  fluids  which  can  be 
claimed  as  in  the  time  period  between  0  to  0.4  s.  The  convection  current  plays  an  important  role  even  though  the  bubbles 
raise  and  then  will  collapse  at  a  shorter  distance.  From  0.4  s  onwards,  the  bubble  growth  takes  place  repeatedly  and  the 
bubbles  move  from  the  surface  freely  and  travel  a  longer  distance.  Since  A1203  Nanofluid  is  having  higher  thermal 
conductivity  compared  to  0.1  and  0.5  VOF,  it  shows  better  HTC  compared  to  Ti02  from  time  period  1.2  s  onwards  there  is 
a  variation  in  the  HTC  for  TiC>2  rather  than  A  ITT,  due  to  the  bubble  density  variation;  the  bubble  raise  will  be  rapid  for  the 
Ti02  Nanofluid  Kshirsagar  et  al.  (2014)  [32]. 

The  following  table  gives  the  initial  and  the  next  preceding  bubble  departure  time.  The  bubble  waiting  time  which 
is  the  difference  between  the  initial  or  first  and  the  second  bubble,  which  allow  us  to  observe  the  time  gap  taken  to  form  the 
next  bubble.  Later  bubble  frequency  which  is  the  inverse  of  bubble  wait  time  is  calculated.  The  results  are  presented  in 
Table  5. 


Table  6:  Bubble  Waiting  Time  (tw)  and  Bubble  Frequency  (f) 


SI. 

No. 

Fluid 

Departure 
Time  (s) 

Second 

Bubble 

Bubble  Wait 
time(tw) 

Bubble  Frequency 
(f)f=l/tw 

1 

W  ater 

3.11 

3.37 

0.26 

3.8461 

2 

Al203-Water 

VOF 

0.05 

1.92 

2.13 

0.21 

4.7619 

0.1 

1.87 

2.04 

0.17 

5.8823 

0.15 

1.78 

1.93 

0.15 

6.6666 

3 

Ti02- Water 

VOF 

0.05 

1.97 

2.16 

0.19 

5.2631 

0.1 

1.62 

1.77 

0.15 

6.6666 

0.15 

1.61 

1.74 

0.13 

7.6923 

As  it  can  be  seen  from  the  above  values,  can  be  concluded  that  higher  the  bubble  frequency,  lesser  will  the  bubble 
waiting  time  which  further  implies  better  TSHF  and  HTC. 

CONCLUSIONS 

To  analyse  the  bubble  nucleation,  growth  and  detachment  from  the  plane,  smooth,  heated  surface  is  carried  out 
using  CFD  the  boiling  phenomena  has  been  simulated  by  an  Euler  description  of  two-phase  flow  (volume  of  fraction 
model).  The  analysed  bubble  profile  show  a  decrease  in  detachment  of  the  bubble  size  and  increase  in  frequency  they  will 
compact  with  each  other  for  the  thermal  energy  from  the  heated  surface.  CFD  results  confirm  the  potential  of  Euler 
two-phase  flow  for  the  boiling  simulation  which  is  useful  to  industrial  applications.  The  studies  concluded  that  by 
increasing  the  volume  of  fraction  for  A1203  and  Ti02  nanoparticles  bubble  frequency  increases  which  can  be  confirmed  by 
the  obtained  HTC  values.  The  studies  reviles  that  after  some  time  of  intervals  Ti02  Nanoparticles  gives  us  the  better  heat 
transfer  rate  compared  to  the  A1203,  due  to  the  increased  in  the  bubble  density  Ti02  VOF  Nanofluid  gives  us  the  better 
results  which  can  be  confirmed  by  the  thermal  physical  properties. 
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